Color-magnetic flux tubes in quark matter cores of neutron stars 

Mark G. Alford* 

Physics Department, Washington University, St. Louis, MO 63130-4899, USA 

Armen Sedrakian"'' 
Institute for Theoretical Physics, J. W. Goethe- University , 
D-6O438 Frankfurt am Mam, Germany 
(Dated: 10 Feb 2010) 

Abstract 

We argue that if color-superconducting quark matter exists in the core of a neutron star, it 
may contain a high density of flux tubes, carrying flux that is mostly color-magnetic, with a small 
admixture of ordinary magnetic flux. We focus on the two-flavor color-superconducting ("2SC") 
phase, and assume that the flux tubes are energetically stable, although this has not yet been 
demonstrated. The density of flux tubes depends on the nature of the transition to the color- 
superconducting phase, and could be within an order of magnitude of the density of magnetic flux 
tubes that would be found if the core were superconducting nuclear matter. We calculate the cross- 
section for Aharonov-Bohm scattering of gapless fermions off the flux tubes, and the associated 
collision time and frictional force on a moving flux tube. We discuss the other forces on the flux 
tube, and find that if we take in to account only the forces that arise within the 2SC core region 
then the timescale for expulsion of the color flux tubes from the 2SC core is of order lO^'' years. 
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I. INTRODUCTION 



It has long been conjectured that neutron stars might contain cores of quark matter, 
and one of the challenges facing nuclear astrophysics is to find signatures by which the 
presence of such matter could be inferred from observations of the behavior of neutron 
stars. This requires us to develop a good understanding of the differences between the 
properties of nuclear matter and quark matter, taking in to account the effects of magnetic 
fields, which are known to be present in neutron stars. In this paper we study quark matter 
in magnetic fields B < 10^^ Gauss, which are astrophysically plausible and high enough to 
affect transport (see for example [1]) but not so large as to modify the phase structure of 
the material [2-4]. 

Nuclear matter at high densities and low temperatures is expected to be a type-II elec- 
trical superconductor, with the magnetic field distributed in an Abrikosov lattice of flux 
tubes [5]. In this paper we investigate the possibility that quark matter in the two- flavor 
color superconducting phase ("2SC") [6] could be a type-II superconductor with respect to 
the color gauge fields [7, 8], with color flux tubes that scatter electrons, muons, and un- 
gapped quarks via the Aharonov-Bohm effect. These tubes are not topologically stable, and 
their energetic stability has not yet been determined; in this paper we investigate the role 
they might play in transport, and their expulsion time, if they turn out to be stable or to 
have a lifetime that is sufficiently long. As we explain below, the tubes carry flux that is 
mostly color-magnetic (hence they can reasonably be called "color-magnetic flux tubes") 
with a small admixture of ordinary magnetic flux. We will argue that the density of color- 
magnetic flux tubes could be high, perhaps only about an order of magnitude less than that 
of ordinary flux tubes in superconducting nuclear matter. Color magnetic flux tubes may 
appear in other color superconducting phases, such as the color-fiavor-locked (CFL) phase, 
but the CFL phase has no gapless charged excitations, so in this paper we focus on the 2SC 
phase. We calculate the Aharonov-Bohm interaction between the flux tubes and unpaired 
quarks or electrons/muons. We calculate the associated damping time and the forces on the 
flux tubes. We defer the calculation of other contributions to relaxation and transport in 
the 2SC phase, such as scattering of the unpaired quarks and electrons off each other, to 
future work. 

The behavior of quark matter phases in magnetic fields is complicated by the intertwined 
breaking of the strong interaction SU (3) "color" gauge symmetry and the electromagnetic 
U{1)q gauge symmetry. In the 2SC phase, a condensate of Cooper pairs of up (u) and down 
(d) quarks leads to the gauge symmetry breaking pattern SU{3)^U{1)q — j- SU{2)rg<S)U{l)Q 
[9, 10]. The unbroken SU{2)rg symmetry ensures confinement of particles that carry net red 
or green color, with a confinement scale around lOMeV [11]. The unbroken U{1)q gauge 
symmetry is a linear combination of the original electromagnetic and color symmetries, called 
"rotated electromagnetism" . The associated gauge field, the "Q photon" , is a combination of 
the original photon and one of the gluons. It is massless and propagates freely in 2SC quark 
matter. The orthogonal combination X is a broken gauge generator, and the associated 
magnetic field has a finite penetration depth. 

The situation is closely analogous to the Higgs mechanism in the standard model, where 
one linear combination of the hypercharge and W3 gauge bosons remains massless (the 
photon), while the orthogonal combination becomes massive (the Z^). The X flux tubes 
are therefore analogous to "Z-strings" [12] which have been found to be stable only in a 
small region of the standard model parameter space [13], although the stable region may 
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be enlarged when bound states are taken in to account [14]. There are differences between 
the 2SC phase of QCD and the Higgs phase of the standard modeh the gluon mass is 
proportional to the quark chemical potential, not the superconducting order parameter [6]; 
the non-Abelian gauge group is SU (3) rather than SU (2) and is only partly broken, leaving 
an unbroken confining SU(2) as well as an unbroken f/(l) in the low temperature phase. 
This means that a separate stability calculation will be needed for the 2SC case. 

Because electromagnetism is much more weakly coupled than the strong interaction, the 
massless Q gauge field is almost identical to the photon, with a small admixture of a color 
gauge boson. Conversely, the broken X gauge field is almost identical to one of the gluons, 
with a small admixture of the photon [10]. Thus the X flux tubes can be described as 
"color-magnetic flux tubes". However, because they contain a small admixture of ordinary 
magnetic flux, they interact with electrons/muons as well as with unpaired (blue) quarks. 
In summary, the 2SC phase is not a superfiuid, but it is_a superconductor with respect to 
the X gauge fields, and a conductor with respect to the Q gauge fields, with current mainly 
being carried by the gapless electrons and blue quarks (one of which is neutral, the other 
has charge +1). Strange quarks and muons, if present, will have a lower Fermi momentum 
because of their higher mass, and hence less phase space near their Fermi surface. Thus 
their contribution to the processes discussed in this paper will be subleading, and we ignore 
it. 

The picture given above is valid below the critical temperature for 2SC pairing and above 
an unknown critical temperature Tisc at which there will be a transition to a phase in which 
there is self-pairing of the blue up and down quarks. Such pairing would break the U{1)q 

symmetry, so there could be both Q and X flux tubes. Models of the strong interaction 
between quarks do not give us much idea of the value of Tisc- They agree that, because 
the strong attraction is much weaker in the single-color channel, Tisc will be many orders 
of magnitude lower than the critical temperature for 2SC pairing, perhaps as low as 1 eV 
(10'* K) [9, 15, 16]. In this paper we will be concerned with temepratures above Tisc, where 
the Q gauge symmetry remains unbroken. 

Depending on the ratio of the X-fiux penetration depth to the coherence length of the 
condensate, the 2SC phase may be type-I or type-II with respect to the X magnetic field [7]. 
In this paper we will be concerned with the possibility of type-II behavior, and the presence 
of flux tubes containing X-fiux in the 2SC quark matter core of a compact star. Even if the 
average magnetic field strength in the core is below the lower critical field, such flux tubes 
may end up "frozen in" if the quark matter had cooled in to the 2SC state in the presence of 
the magnetic field. The magnetic field would then be resolved in to a Q part, which would 
pass freely through the 2SC quark matter, and an X part, which would become trapped in 
flux tubes (Sec. III). 

The paper is structured as follows. In Sec. II we calculate the Ginzburg-Landau parameter 
for 2SC quark matter, and conclude that it is a type-II superconductor with respect to 
the broken X generator as long as the pairing gap A is large enough. We estimate that 
A > /ig/16 will suffice, which for typical quark chemical potentials fig ~ 400 MeV requires 
A > 25 MeV. In Sec. Ill we discuss the nucleation scenario by which the flux tubes can 
occur in the 2SC superconductor, even when the magnetic field intensities are below the 
lower critical field. We estimate the density of such flux tubes in the hypothetical 2SC quark 
matter core of a neutron star. In Sec. IV we calculate the Aharonov-Bohm scattering cross 
section for electrons or unpaired quarks interacting with color magnetic flux tubes. Sec. V 
is devoted to the computation of relaxation time of massless electrons and unpaired blue 
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quarks interacting with flux tubes via Aharonov-Bohm cross-section. In Sec. VI we estimate 
the timescale for expulsion of the flux tubes from the 2SC core, taking in to account the 
forces on the color-magnetic flux tubes in the 2SC core and at its boundary, but neglecting 
any forces on the magnetic flux lines outside the core. We summarize our results in Sec. VII. 

In our calculations we use "Heaviside-Lorentz" natural units with h = c = ks = = 1, 
where is the Boltzmann constant and eo is the vacuum permittivity; the electric charge 
e is related to the fine structure constant by a = e^/(47r). 



II. TYPE-II COLOR SUPERCONDUCTIVITY IN QUARK MATTER 

A superconductor is of type II if it obeys the condition 

where k is the Ginzburg-Landau (GL) parameter, A is the penetration depth, and ^ is the 
coherence length for the superconductor. In a system of relativistic fermions with chemical 
potential and pairing gap A, we expect ^ oc 1/A, A oc {gfi)~^, so k oc A/{gfi). (In the case 
of 2SC quark matter the relevant broken gauge symmetry is the "X" which is mostly color, 
so the coupling g is approximately the strong coupling constant.) We therefore expect that 
2SC quark matter will be a type-II color superconductor if the gap is sufficiently large. 

To make a more accurate determination we follow the approach of Bailin and Love [17] 
and lida and Baym [7]. We start with the effective free energy density (Ginzburg-Landau 
theory) for a relativistic BCS superconductor (Ref. [17], (3.12)) 

= + aij*^ + + 7(VV^* - 2ieAij*){Vij + 2ieA^) + ^{B - ^i^Hf . (2) 

2 2/io 

(We have followed Ref. [17] in writing the magnetic field free energy in SI units; in natural 
units /io = 1.) Here ip is the gap parameter; for negative a the free energy has a minimum 
at l^/'p = ipQ, with penetration depth A, and coherence length ^ given by 

^0 = 4' A^= , e = -^-. (3) 

where gpair is the charge of the Cooper pair. The GL parameter k is then given by 

2 1/3 

The coefficients in the Ginzburg-Landau functional are [17] 

tgl 

^ = ^lHw' (5) 
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where tgl = {T — Tc)/Tc. The fermions have Fermi momentum pp, so the density of states 
near the Fermi surface is = Nppfi/n^ ~ Nfi'^/n^. The parameter is a degeneracy factor 
that is 1 for a single-species system, and 2 for the 2SC phase (see Ref. [17], Eq. (4.63)). The 
Ginzburg-Landau theory is most rehable for temperatures close to Tc, however we will use 
it at T ^ Tf,. The low-temperature gap parameter A is related to the critical temperature 
by Tc = {e'^/7r)A: note that A then differs from ipQ by a factor of about 1.7. Expressing the 
coefficients in terms of A, 

32.74 A 

(6) 



gpairVA^ 

We can check this result by noting that for a relativistic electronic superconductor, = 1 
and gpair = 2e where a = e^/(47r) ^ 1/137. Substituting these values into (6) we find 
K = 54.043A//i = 95.325Tc//i, in agreement with Ref. [17], (3.24). 

In 2SC quark matter, the degeneracy factor is = 2 and the charge of the Cooper pair 
is the X charge of the 2SC condensate. From Eqs. (29) and (31) of Sec. IV we find 

v3cos(p V3 

where the mixing angle (p is defined in Eq. (29). We estimate the strong coupling constant 
g by assuming that = g'^/{47r) ^ 1, so g ^ 3.5. Substituting these values into (6) we find 

«2sc ~ 11— • (8) 

We conclude, using (1), that 2SC quark matter will be of type II if the pairing gap is 
sufficiently large, A > fiq/16. In quark matter we expect fig ~ 400 MeV, so this only 
requires the 2SC pairing gap to be greater than about 25 MeV, which is well within typical 
estimates [18, 19]. Our general conclusion agrees with that of Refs. [7, 8] who also noted 
that a sufficiently large 2SC pairing gap yields a type-II superconductor. Our specific result 
(8) differs from Eq. (112) of Ref. [7] by a factor of a/2, but given the uncertainty in the 
strong coupling constant g this numerical discrepancy does not affect our conclusion. 



III. COLOR-MAGNETIC FLUX TUBES IN THE 2SC PHASE 



A. The nucleation and density of flux tubes 

When the quark matter core of the star cools below a critical value, a 2SC condensate 
forms. We expect that this happens before the nuclear mantle becomes superconducting 
because the gap parameter for quark matter is expected to be an order of magnitude larger 
than that for proton pairing [6, 20-22]. The electromagnetic field is then resolved in to a 
Q component and an X component. The 2SC core is not a superconductor with respect 
to Q, so the Q component is undisturbed [10] (on this we disagree with Ref. [23], which 
we believe imposes an incorrect boundary condition on the gluon field). However, the core 
is a superconductor with respect to the X component, and we have argued above that it 
may well be a type-II superconductor. The lower critical field for the X-superconductivity 
is very high, Hd ~ 10^^ Gauss [7], and typical neutron star magnetic fields are expected to 
be lower than this, but, as we now argue (see also [10] and footnote [8] of Ref. [7]), it is still 
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quite possible for the X-flux to form flux tubes threading the quark matter core. The only 
way the X-flux could be expelled from the core is if the transition from hot quark matter to 
2SC happens smoothly from the center of the star outwards. However, it seems more likely 
that the transition to 2SC matter will proceed by nucleation of 2SC regions ("bubbles") in 
the quark matter, which then grow and coalesce. The X-fiux will be expelled from the 2SC 
bubbles, but will then be trapped in the non-superconducting regions between the bubbles. 
As the bubbles grow, these regions become smaller, concentrating the flux there until the 
local field strength rises above Hd, at which point the bubbles stop growing. At this stage, 
the core consists of 2SC quark matter with channels of non-superconducting quark matter 
running though it, carrying the X-fiux. If the 2SC phase is a type-II superconductor then 
these channels are unstable and will fragment into flux tubes, each carrying a single quantum 
of X-fiux, with a short-range repulsion between the flux tubes. The fact that the average 
field strength was below the lower critical field for a sphere of 2SC matter in a uniform 
magnetic field will now manifest itself as an outwardly-directed boundary force on the flux 
tubes at the point where they meet the edge of the 2SC core. We will study this in Sec. VI. 

Because the 2SC phase is a conductor with respect to Q charge, it supports eddy currents 
which make it very difficult for the Q magnetic field in the 2SC core to change. The timescale 
for expulsion of the Q magnetic field is estimated to be longer than the age of the universe 
[10]. Thus we are justified in treating the Q magnetic field as a fixed background. 

If we assume that all the X-fiux is trapped in the manner described above, then the 
density of flux tubes is just Bx, the density of magnetic X-fiux, divided by $x, the X-fiux 
of a single flux tube. Bx is obtained by projecting out the X-component of the original 
electromagnetic flux B (see (27)), so Bx = Bsiinp. The flux quantum is 



<l>x = -\—^ (9) 

'Zpair V 

where gpair is the AT-charge of the 2SC condensate (see (7) and (29)). We can relate it to 
the flux quantum = ir/e ^ 10.37 for an ordinary superconductor where the charge of the 
condensate is 2e, 

2e 

$x = $0 = 6sin(y?)$o . (10) 

Q'pair 

We conclude that 

This is the upper limit on the flux tube density in 2SC matter. Interestingly, as anticipated in 
Ref. [24], it only differs by a factor of 1/6 from the density of electromagnetic flux tubes that 
would result if the core were an electromagnetic superconductor due to electron or proton 
pairing. Projection on to the X component reduces the magnetic flux by a factor sin if, but 
because the X flelds are strongly coupled their flux quantum is smaller by a similar factor, 
so the flux tube density ends up being independent of the mixing angle. The actual density 
will depend on details of how the transition to 2SC matter was completed. For an internal 
fleld B = 10^^ Gauss (2MeV^), the maximum flux tube density is = 8.1 x lO^^cm"^. 



B. Properties of the flux tube 

The thickness of the flux tubes is given by the penetration depth for magnetic X-flux in 
the 2SC phase. This follows from equations (4) to (7). Assuming pp — fiq for relativistic 
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quarks, 

The energy per unit length (tension) of the flux tube is given by \E In kx where E is the 
energy per unit length of the magnetic flux if it were uniformly spread over an circle of radius 
A (Ref. [25], Sec. (5.1.2)), and In^x is a factor of order 1. In Heaviside-Lorentz natural units 
E = By 2, where B = ^^/(vrA^), so 

(compare Ref. [7], Eq. (107); see also Ref. [24]). To estimate the tension we work to lowest 
order in a and use (12), (9), and (34). In the low temperature limit we find 

u2 

ex = !^\nKx. (14) 

Assuming that in 2SC quark matter fig is in the 350 to 500 MeV range, and that the loga- 
rithmic factor is of order 1, we conclude that the tension will be of order 60 to 130 MeV/fm. 



IV. AHARONOV-BOHM SCATTERING BY FLUX TUBES 

The Aharonov-Bohm effect provides a remarkably strong interaction between a charged 
particle and a flux tube containing magnetic flux. For the simple case of a single U{1) gauge 
group (electromagnetism) , the differential cross-section per unit length is (see, for example, 
Ref. [26]) 

da _ sin^(7r/3) 

where 

/3 = ^ , (16) 

where qp is the charge of the scattering particle. For a flux tube that arises as a topological 
soliton in an Abelian Higgs model, Qc is the charge of the condensate field whose winding by 
a phase of 27r characterizes the flux tube; k is the momentum in the plane perpendicular to 
the string, and i) is the scattering angle. Aharonov-Bohm scattering has several important 
features: 

• The cross-section vanishes if /3 is an integer, but is otherwise non-zero. 

• The cross section is independent of the thickness of the flux tube: the scattering is not 
suppressed in the limit where the symmetry breaking energy scale goes to infinity, and 
the flux tube thickness goes to zero. 

• The cross section diverges both at low energy and for forward scattering. 

It is therefore of great interest to determine the values of /3 for scattering of the fermions 
that are ungapped in the 2SC phase off a flux tube containing magnetic flux associated with 
the broken gauge symmetry. 
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A. The gauge groups and charges 



1. The light fermions 

In the 2SC phase we will focus on the U{1) x U{1) gauge group consisting of electromag- 
netism and the part of the color gauge symmetry that mixes with electromagnetism. The 
relevant particles are the quarks and the electron: 

ip = {ru, gd, rd, gu, bu, bd, e~), (17) 

where "rw" means the red up quark, etc. "e~" is the electron. Muons would have the 
same interaction as the electron, so we do not include them separately. In this basis, the 
generators of the two U{1) gauge groups are just the diagonal matrices of their electric and 
color charges, 

= diag(+|,-|,-|,+|, +1,-1,-1), 

= i73diag(l, 1,1, 1,-2, -2,0). ^ ^ 

The normalization of Q"^ is fixed by the conventional electric charges of the particles. For 
T we have used the conventional normalization for generators of the SU{3) color gauge 
group [27]. The kinetic term in the lagrangian of the fermions is ip'y'^D^ip, where the covari- 
ant derivative of the fermion fields is 

D^^ = d^^ - leA^Q^^ - igA^T^ . (19) 

The electromagnetic gauge coupling is e, and the QCD gauge coupling is g. The photon 
gauge field is A'^, and the gluon gauge field is A"^. With the normalization of (18), a = 
6^/4:71 = 1/137, and a, = g^/iTi ~ 1. 



2. The 2SC condensate 



The 2SC condensate is a diquark condensate, 

= (tpiC'y^ipj) 



(20) 



where the indices i and j live in the color- flavor space of (17). The condensate only involves 
the red and green up and down quarks, so its color-flavor structure is 
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(21) 



From (18) we can see how (pij, considered as a 7 x 7 matrix in the color-flavor space of (17), 
transforms under an infinitesimal electromagnetism or color rotation. Each of the quarks in 
the diquark feels its own color-flavor phase, so each index i and j is separately transformed: 



Q'f' ■(j) + (j)-Q^, 
■ + ■ T'^, 



(22) 
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where the "■" on the right hand side signifies ordinary matrix multiphcation of the two 
7x7 matrices, and we have used the fact that and are both diagonal, and hence 
symmetric. 

The lagrangian of the 2SC condensate {i.e. the G-L theory) contains the kinetic term 
{D^(j))*D^^(j), where the covariant derivative is 

D^cP = 9^0 - 2eA«Q^0 - igApU ■ (23) 

This determines the couphng of the 2SC condensate to the gauge fields. 

B. The broken/unbroken basis 

When the 2SC condensate forms, one linear combination of Q and T is spontaneously 
broken: we will call it "X" . The other remains unbroken: we will call it "Q" , 

X = -r]2Q + T . 

We determine r]i by requiring that the 2SC condensate be invariant under Q gauge trans- 
formations, 

QU = , (25) 

which implies that 

"71 ■ ^''^ 

It is natural to work in the {Q, X) basis rather than the (Q, T) basis, so we define new 
"rotated" gauge fields 

A^ = simpA'^ + cos(pA^, ^ ' 

where the mixing angle ip is analogous to the Weinberg angle in the standard model which 
parametrizes the mixing of the hypercharge and gauge bosons to yield the photon 
(analogous to A^ here) and the Z (analogous to A-^ here). It is important that the mixing 
of the gauge fields is expressed in terms of an angle, so it maintains their normalization, so 
the gauge field kinetic terms for A'^ and A^ are still conventionally normalized. In the case 
of the generators, which we defined in (24), the overall normalization is not important, since 
it is absorbed in the new gauge couplings. 

In the new basis, the covariant derivative of the fermions is 

D^iP = d^^ - le^^^A'^^Q'f'^ - ie^^MJX^V' • (28) 

We will determine the new gauge couplings e*^'^-' and e^^\ and the mixing parameters 772 and 
by requiring that (28) be equivalent to (19) for all gauge field configurations. 
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C. X-charges of the particles and condensate 



Flux tubes will contain magnetic X-fiux, so to determine the Aharonov-Bohm scattering 
parameter P for each particle, we need to find the X-charge of each particle, corresponding 
to Qp in (16). This follows straightforwardly from (28). We will also need to know the 
X-charge of the 2SC condensate, corresponding to Qc in (16). 

Requiring that (28) be equivalent to (19) for all gauge field configurations, and using 
(26), we find 

-V^ tan^ ip 

(29) 



cos 



g2 

m = ^ — = — V^tan^o9 



,(Q) 



y/Seg 

ecosip 



<x) 

' = — , = g cos 09 . 

There is a new "rotated" electromagnetism, with coupling e'-'^-', which is slightly smaller 
than the usual electromagnetic gauge coupling. The charges of the fermions under this 
gauge group are 

g = diag(i,-i,-i,i, 1,0,-1) . (30) 

This agrees with the well-known results for the 2SC phase [6]. 

The action of the X-charge matrix on the 2SC condensate determines the X-charge qc of 
the condensate, in units of e^^^; from (21), (24), and (29), 

X0 + 0X = q^(f), 



where Qc = (l + 



(31) 



-\/3^ 3(7^/ \/3cos^Lp 



=2 



The X-charge matrix of the fermions is 

1 

73 

2 



X = ^=diag(| + 2 tan^ v?, | — tan^y?. 



— tan^ (p, 1 + 2 tan^ ip, 
-1 + 2 tan^ p, —I — tan^ (p, 

2 



(32) 



-3 tan Lp). 



Dividing by Qc (31) we find the Aharonov-Bohm /3-factors of the fermions, in the basis 
defined by (17), 

= diag(^ ^ + ^ sin^ \ ~\ ^^^^ 
1 3 . o 13 



2,„ - , -„-„2,„ (33) 



sm c/9, — — sm p 

2 2 ^22 ^ 



-l + 3sin^V9, —1, — 3sin^</?j 
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Expanding in powers of (since e ^ (7), we find 



sin^(y.) ^ — (34) 

6a. 



so to lowest order in a, 

^'/' = diag(^ + 



1 a 1 a 



2a, 2 2a, 
1 a 1 a 



2 2a, 2 2a, 



(35) 



a ^ a 
+ — , -1, 

a, a. 



We conclude that the gapped quarks have (3 close to |, which means that they have near- 
maximal Aharonov-Bohm interactions with an X-fiux tube. Among the lighter (and hence 
more phenomenologically relevant) fermions, the Q-neutral bd has zero Aharonov-Bohm 
interaction with the flux tube, while the bu and electron have the same Aharonov-Bohm 
factor 

sin(7r^''") = sin(7r;g") ^ -n— . (36) 

a, 

V. RELAXATION VIA SCATTERING OFF FLUX TUBES 
A. Relaxation time calculation 

In this section we compute the characteristic timescale for a perturbation from equilibrium 
to relax away due to scattering of the fermions off the color magnetic flux tubes. This 
relaxation time is a measure of the mean free time between collisions of the fermions with 
the flux tubes, so we will also refer to it as a collision time. Our calculation applies equally 
to electrons and the unpaired component of blue colored quarks in the 2SC phase, the key 
difference being /3 factors in the cross- sect ion. The Boltzmann kinetic equation for blue- 
quark/electron distribution function f{p,t) is 



dfip,t) 2ixN„ fdy 



W{p;p')f{p',t) [l-/(p,t)] 



dt V J (27r)^ 

- Wip';p)fip,t)[l - fip',t)]jS{e{p) - e{p')), (37) 

where A^^, is the number of flux tubes, V is the volume, and W{p';p) is the transition 
probability between the states described by momenta p and p' . Time-reversal symmetry 
implies W{p]p') = W{p']p). In equilibrium the fermion distribution function is given by 
the Fermi-Dirac distribution function 

Up) = T- rj^ 777^ (38) 

1 + exp[(p-/ii)/T] 

where T is the temperature and /Xj is the chemical potential of blue quarks {i = b) and 
electrons {i = e). To solve the Boltzmann equation we shall apply the variational method, 
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where the perturbations from equihbrium are described by variational trial functions whose 
functional form is dictated by the form of applied perturbation [28, 29]. The resulting 
transport coefficients are lower bounds on their exact values. The number of adjustable 
trial functions, which are used to maximize the entropy production via scattering, could 
be large. In the following we shall use one linear function 0, in which case there is no 
need for variation, since the variational parameter cancels out. It should be kept in mind 
that the resulting transport coefficients are still lower bounds on their exact values. For 
small perturbations from equilibrium the Boltzmann equation can be linearized by writing 
f{p,t) = fo{p) + Sf{p,t), where the (small) perturbation from the Fermi-Dirac form (38) is 

Sf{p,t) = -^<P{p,t), (39) 

where (plp, t) is the trial function. The linearized Boltzmann equation then reads 

-^/o(p) 11 - /«(P)1 = ^ /(gs Wp.')-*(p'.')1 

X W(p,p')h(p') |1 - /„(p)l i(e(p) - e(p')). (40) 

To obtain this form of the kinetic equation we used the detailed balance conditions 
Up') [1 - Up)] - Up) [l - Up')] = O, and dUp)/de{p) = dUp')/de{p'). It is conve- 
nient to work with the Laplace transformed trial function 

0(p,t)= / rfse"^V(P,s). (41) 



Upon Laplace transforming Eq. (40) we find 

s<P{p,s)Up)[i-Up)] = ^ J0^s[<l^(p^^)-<f>ip'^^)] 

X W{p; p')Up') [1 - Up)] mP) - e{p')). (42) 

To define a characteristic relaxation rate we assume that the trial function can be written 

as 

<f){p,s) = <f){p)6{s-so), (43) 

in which case Eq. (42) becomes 

soHp)Up)[i-Up)] = ^ j^Mp)-^^p')] 

X w{p-p')Up') [1 - Up)] mp) - e{p')), (44) 

where the perturbation functions are now independent of s. We can identify Sq with the 
relaxation rate (z.e.the inverse of the relaxation time) by comparing the computed kinetic 
coefficients with standard expressions for transport coefficients, e.g. , the electrical conduc- 
tivity with the Drude formula. 

To formulate the variational principle [30] we write Eq. (44) in the compact form 

X{p) = J [<P{p) - <jy{p')] P{p, p')d'p', (45) 
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where X{p) stands for the left-hand side of Eq. (44); the scattering operator P(p, p') is easily 
read-off from the kernel on the right-hand side of Eq. (44). Since the factor /o(p) [1 — /o(p)] 
and the transition probability are positive definite, the operator P{p,p') is positive definite. 
Further it is linear and self-adjoint (symmetric). Following ref. [30] we define an inner 
product 

= j 0(p)^(p)rfp, (46) 

in terms of which 

(0,P^) ^IJdpJ dp'[<P{p)-<P{p')]P{p,p')mp)-i;ip')]. (47) 

The variational principle states that the expression 

(0,X) = (0,P0) (48) 

attains its maximum for the exact value 0ex; which satisfies Eq. (45); for any other trial 
function (j), that satisfies Eq. (48), {(l),P(j)) < (0cx,P0cx)- Exphcitly, Eq. (48) reads 

I (0s*(p.-')V«(P)ii-/.(p)i - ^ /g,iwp..)-*(p',.)f 

W{p;p')Mp') |1 - /„(p)l *(£(p) - £(p')). («) 

It is also straightforward to check that the variation of Eq. (49) leads us back to "equation 
of motion" (44). From Eq. (49) we obtain the variational relaxation rate 

(50) 

where 

1^ = 1 ^.m'foiP) [1 - /o(p)] . (51) 

The exact relaxation rate s > Sq. We specify the form of the trial function appropriate to 
the problem at hand, which is the relaxation of uniform blue-quark/electron velocity v on 
a flux tube, as 

^{p)=p-vC{p'), (52) 

where C{p^) is the scalar part of the trial function. In the following we will adopt the simple 
choice C(p^) = 1. The differential transition probability can be obtained from the Aharonov- 
Bohm scattering cross-section, Eq. (15), and is given by (for details see Appendix A) 

,53) 

where the initial and final state momenta and energies, p and e, are unprimed and primed 
respectively (in this section we use for the scattering angle, as opposed to d in Sec. IV). 
Here we have used the cylindrical coordinates coaxial with the flux tube to write d^p = 
p±dp±d(f)dpz. Combining Eqs. (50) and (53) and carrying out the phase space integrals 
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(the details are given in Appendix B) we obtain, to lowest order in the low-temperature 
expansion, 

where ppi is the blue-quark/electron Fermi momentum, is the density of flux tubes. 
Eq. (51) with the trial function (52) can be computed in the low-temperature limit by 
approximating df{p)/de{p) ~ 6{e{p) — fj,i) to obtain V = pp^v'^T/Gn'^, where v is the fermion 
fluid velocity (52). The relaxation rate for particles of species i scattering off flux tubes of 
area density n„ is then given by 



-1 IT'v . 2 



if 



Sq = — sin (7r/3i) . (55) 



PFi 



It easy to understand the final result (55). It is of the standard form for classical gases 
= cna, where c = 1 is the speed of the particles, n = n^, is the density of scattering 
centers, and a oc sm"^ {n 13) /pf is the cross section for Aharonov-Bohm scattering. Eq. (55) 
is relevant for thermal relaxation of the gapless fermion species in the 2SC phase. One of 
these, the blue down quark, has no A-B interaction with the flux tubes (/3 = 0). The other 
two, the electron and blue up quark, have identical A-B factors (36) although their Fermi 
momenta are different. 



B. Comparison with Coulomb scattering 

To find out whether scattering off flux tubes is likely to be an important source of re- 
laxation, and hence a significant contributor to transport properties, it is useful to compare 
Eq. (55) with the collision time for screened Coulomb scattering via exchange of Q photons. 

The 2SC phase is a Q-conductor, with two species of gapless charged fermions^ the bu 
quarks (with Q-charge +1 and chemical potential ~ fi) and the electrons (with Q-charge 
-|-1 and chemical potential ~ fie)- There may also be muons, but their Fermi momentum 
will be much smaller. As mentioned in the introduction, the red and green quarks will be 
confined to bound states whose mass is expected to be of order lOMeV [11], so they play 
no role in transport at neutron star temperatures. Since fi > fie, the bu quarks are more 
numerous than the electrons and have a larger phase space near their Fermi surface, so they 
will make the largest contribution to the collision time. 

The Coulomb collision time depends on the in-medium photon spectrum, which will 
be affected by Debye screening and Landau damping arising from the presence of gapless 
charged excitations, dominantly the bu quarks because of their larger phase space. A simple 
estimate can be obtained by assuming that the dispersion relation is dominated by a plasmon 
pole. The plasma frequency Up is given by 

-v = — = ^.<^N (56) 

~ 2 ~ 

where a = e^*^^ /(4vr) is the fine structure constant for the "rotated" Q electromagnetism 
(29). The colhsion frequency is given by (see Eq. (10), (12), (18) of [31]), 

_i 803)^ _ 6C(3) 
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where ({3) = 1.202. This result is vahd for T ^ Up, which is the relevant regime for neutron 
stars since fig is in the 400 MeV range. Eq. (57) is analogous to Ref. [32] 's Eq. (62) for the 
thermal conduction timescale, with electromagnetic interactions (so their as is replaced by 
(5) and a different number of quark species. The quark-quark Coulomb collision frequency 
(57) is proportional to temperature T whereas the particle-flux-tube collision frequency is 
independent of temperature. We can therefore define a temperature Tf below which flux 
tubes dominate the relaxation of deviations from thermal equilibrium. From (57) and (55) 
we find 

To make a numerical estimate we assume that the 2SC core contains the maximum flux tube 
density given by Eq. (11), and that as ~ 1. Using (36) and the fact that a ^ a, we find 

We conclude that for reasonable values of the magnetic field, only at very low temperatures 
is Aharonov-Bohm scattering off flux tubes likely to be an important source of thermal 
relaxation. However, it is important to note that the thermal relaxation timescale is not the 
only one that is relevant to transport. There is also the viscous relaxation rate (Ref. [32], 
Eq. (51)) and the momentum relaxation rate (Ref. [32], Eq. (32)) both of which have a much 
stronger (oc T^/^) suppression at low temperatures. We defer a full discussion of transport 
in the 2SC phase to later work. 



VI. FORCES ON THE FLUX TUBES 

We argued in Sec. Ill that even if the magnetic field in the core of the star is below 
the lower critical field, color magnetic flux tubes will still be produced in the transition 
to the 2SC phase. In this section we study the forces on those flux tubes, and estimate 
the timescale for their expulsion from the 2SC core. For this initial estimate we take in to 
account only the forces on the flux tubes within the 2SC core, or at its boundary. Depending 
on the nature of the material surrounding the core there may be additional forces, and these 
may modify the expulsion time in a way that would have to be calculated on a case-by-case 
basis. 

The velocity of the flux tube is Vl, the velocity of the normal fluid is v^, and the 
velocity of the 2SC condensate is vs- The forces we consider are mutual friction ("mf"), the 
non-dissipative (lifting) Magnus-Lorentz force ("ML"), the lordanskii force ("lord"), forces 
arising from zero modes ("zm"), and boundary forces ("bf") at the quark-hadronic boundary. 
We assume that local magneto-hydrostatic-gravitational equilibrium is established quickly 
after the transition to the 2SC phase, so there is no additional buoyancy force [33, 34]. We 
note that there may be additional forces due to density-dependence of the 2SC pairing gap 
[35], but we do not include these since there is as yet no reliable estimate of the density 
dependence. The equation of motion of a flux tube then has the form 

'^^~ir = -^mf + -^ML + /lord + /zm + /bf , (60) 
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where my is the effective mass of a flux tube per unit length and / is a force per unit length. 
The boundary forces tend to pull the flux tube in a radial direction, expelling it from the 
2SC core. This is resisted by the combination of the other forces. 

In our calculations we will assume that the flux tubes are straight. A bent flux tube will 
feel an additional restoring force determined by its tension. 



A. The background Q magnetic field 

In our calculations we will neglect the effect of the Q magnetic field Bq that penetrates the 

2SC core. Because of this field, Q-charged particles, including the bu quarks and electrons, 
will feel a Lorentz force. This will have a significant effect on the behavior of the normal fluid 
of quarks and electrons when the cyclotron frequency Uc becomes larger than the inverse of 
the characteristic time for equilibration, which as we argued in Sec. VB, is the quark-quark 
collision time Tgg (57). The dominant component of the fluid is the bu quarks, with Q-charge 
~ e, and Bq ^ B, so 

= — , (61) 

Pf 

and we can neglect the effects of the magnetic field on transport when UcTqq <^ 1, where 

271^ 1 B f B \ /lO^Kx /400MeV\ 

We conclude that only for high magnetic fields (above 10^^ G) or low temperatures (below 
10^ K) might the magnetic field affect thermal relaxation. We defer a discussion of this 
regime to future work. 



B. Mutual friction 

Mutual friction is a frictional force on a flux tube arising from its Aharonov-Bohm in- 
teraction with the normal fluid of gapless particles through which it is moving. Consider 
a vortex moving relative to the normal fluid with velocity u = Vl — Vn- In the relaxation 
time approximation 

^-^ = ^/(^^^°^^'^^' ^^'^ 

where r^j^ is the collision rate between fermions of species i and flux tubes (55). 

We will assume that the blue up quarks dominate the friction. This is reasonable because 
the blue down quarks have no X charge and hence no Aharonov-Bohm interaction with the 
flux tube, and the electron Fermi momentum is smaller than that of the blue quarks. The 
equilibrium Fermi-Dirac thermal distribution of the quarks is 

/o(p, u) = {exp[(£ - /i. + p ■ u)/T] + 1}-' , (64) 

where the p u term is a correction due to the motion of the vortex relative to the thermal bath 
with velocity u. We will compute the force to linear order in u. The leading contribution 
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FIG. 1: A straight flux tube of lengtfi I passing through a 2SC neutron star core of radius R, at 
distance r from the center of the star. There is a boundary force where it reaches the edge of the 
2SC core. 



arises at the first order in velocity 



1 r ^3 



/mf = — / J^PiP ■ U)^^^^^ = 7]U. (65) 

J [2'kY de 

Carrying out the integral and using (55) we obtain the mutual friction drag coefficient 

^ = ^^' ' =n,sin^(7rA) (66) 

Til) 

where rij is the fermion density and is their Aharonov-Bohm factor (35). As one would 
expect, the friction coefficient is independent of the magnetic field (i.e. the density of flux 
tubes). It is proportional to the fermion density, so, as noted above, the hu quark contribu- 
tion will dominate the electron contribution. 



C. Boundary forces 

Next we wish to calculate the force exerted on the flux tube at the point where it reaches 
the interface between the 2SC quark matter core and the nuclear mantle of a neutron star 
(Fig. 1). When the X-magnetic flux tube reaches the edge of the 2SC core it combines with 
the Q magnetic flux in the core to re-constitute the ordinary magnetic fleld from which it 
was originally formed. The form in which the flux continues through the nuclear mantle, 
and hence the boundary energy, may therefore be influenced by the state of the nuclear 
matter. In this analysis we will include only the forces arising from the contribution due to 
the 2SC core itself. We briefly discuss other contributions, but a proper analysis including 
them would have to be done in the context of a speciflc model of the whole neutron star 
and the properties of all regions within it. 
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The outward force per unit length on the flux tube is (see Fig. 1) 



IdE r , _ 

^' = 1^ = ^^"^ • ^^^^ 

Here E should be the total energy of magnetic flux inside and outside the core, but we 
neglect the outside contribution; ex is the energy per length of X flux tubes. Then from 
(14), 

/"^^^T^l^'^^ • (68) 

Taking in to account the energy of the magnetic fleld outside the core will reduce the 
right hand side of (67), and weaken the outward force on the flux tube. We now discuss the 
magnitude of such terms in various cases. 

If the nuclear mantle is a type-II superconductor, the magnetic fleld penetrates the nuclear 
mantle in the form of Abrikosov flux tubes (dashed line in Fig. 1). From (11) we know that 
each X flux tube will spawn 6 Abrikosov flux tubes in the nuclear matter, each of which has 
energy per unit length 



$2 







luKnuc • (69) 



^0 = 7r/e ~ 10.37, so if we assume that the logarithmic factor is of order 1 then for Anuc in 
the 50 to 100 fm range, ^nuc is in the 0.2 to 0.7 MeV/fm range. This means that even when 
multiplied by a factor of 6, ^nuc is small in comparison with the tension of the X flux tube, 
which is greater than 10 MeV/fm (14), so (68) is still a good estimate of the boundary force. 
Of course, in a type-II nuclear mantle there may be other forces, for example if it is also a 
superfluid there may be entanglement of Abrikosov flux tubes with superfluid vortices, but 
we neglect those here because they depend on details of the nuclear mantle. 

If there is no Cooper pairing of the protons then the nuclear matter is a conductor. In this 
case the energy gained from shortening the X flux tube is counteracted by the fleld energy of 
the magnetic fleld it connects to in the nuclear matter mantle. The criterion for the tension 
of the flux tube to dominate is the same as the criterion for the magnetic fleld to be below 
its lower critical value. Since neutron star magnetic flelds are well below the lower critical 
fleld for the 2SC phase, we can assume that the 2SC flux tube tension will dominate and we 
can use (68) again. The only complication is that conducting nuclear matter supports eddy 
currents which will resist any change in the magnetic fleld in the nuclear mantle. This may 
make it much harder to move the X-flux tubes in the 2SC core. Again, we do not attempt 
to include such forces that depend on details of the constitution of the nuclear mantle. 

If the nuclear mantle were a type-I proton superconductor [36-39] then X flux tubes 
in the 2SC core would connect to non-superconducting domains in the nuclear mantle [36, 
39]. In this case we cannot compute the boundary force because the domain structure 
of the type-I proton superconductor is not known; the possible (layered, cylindrical, etc) 
structures in type-I superconductors essentially depend on the history of the nucleation of 
the superconducting phase. 



D. Magnus-Lorentz force 

The Magnus-Lorentz force is a non-dissipative force, directed orthogonally to the flux tube 
velocity, that arises from the superposition of the winding "flow" of the 2SC order parameter 
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around the flux tube and the background flow of the charged superfluid of fermions [40-42] . 
(There is controversy about this in the literature; for example, Jones [33, 43] has suggested 
that this is cancelled by another contribution from ungapped fermions. Pending a definitive 
resolution of this disagreement we will use the standard form of the Magnus-Lorentz force.) 
The Magnus-Lorentz force per unit length on a flux tube is 

/ml = -{jx X n^x) , (70) 

where $x is the X-flux through the flux tube (9), n is a unit vector pointing along the flux 
tube, and jx is the current of X charge seen by the flux tube, arising from the X charge 
density px of the 2SC condensate, moving relative to the flux tube 

jx = pxivs - vl) ■ (71) 

We can write px = where Ug is the density of quarks in the condensate. Since there 

are 4 quark species in the condensate, and at low temperature all fermions are part of the 
condensate, 

/ml = -pivs -VL)xn, 

4;.3 (72) 
P = Px^x = T^ris = —— . 

371 

Note that the charge of the Cooper pairs cancels in this expression. 



E. lordanskii force 

The mutual friction force described above is the force on the flux tube in the longitudinal 
direction (i.e. parallel to its velocity relative to the normal fluid of unpaired quarks), due to 
Aharonov-Bohm scattering of the unpaired quarks. The lordanskii force is the transverse 
component of that same force [42] , 

/lord = D' {vl -VN)y<n. (73) 

The transverse Aharonov-Bohm scattering cross-section for bu quarks off the flux tube is 
(T_L = — A;~^ sin(27r/3''") (Ref. [42], Eq. (64)) and, as in the case of the longitudinal Aharonov- 
Bohm force, one expects the force per unit length to be proportional to the fermion density, 
so we expect D' ^ sin(27r/3''")/i^ ^ a/x^ (see Ref. [42], after Eq. (69)). This rough estimate 
is sufficient to argue that the lordanskii force can be neglected. Basically, the Aharonov- 
Bohm forces are suppressed by powers of a arising from the Aharonov-Bohm factor of the 
bu quarks (36). In the case of the lordanskii (transverse) component, we will see that this 
makes it subleading relative to the Magnus-Lorentz force, which also acts perpendicular to 
the flux tube's velocity. In the case of the longitudinal component, there is no larger force 
parallel to the velocity, so the Aharonov-Bohm force is the dominant contribution to mutual 
friction. 



F. Zero-mode force 

The frictional force on a flux tube due to scattering of zero modes localized inside the 
flux tube off gapless fermions in the bulk [44, 45] has been calculated for proton flux tubes 
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in nuclear matter [33]. At low temperatures, we expect the frictional force on a 2SC flux 
tube to be 

/|| = -^K-M, (74) 

where, generalizing from nonrelativistic protons to relativistic quarks, 

C = Trrig tanh(A/2T) ~ /xj , (75) 
uo ~ AV/i, , (76) 
Tc - /.f . (77) 

Eq. (75) follows from Ref. [33] Eq. (7), and the fact that 2SC pairing gap A is expected to be 
much bigger than typical neutron star temperatures. Eq. (76) follows from Ref. [33] Eq. (1), 
assuming, following Ref. [33], that the typical transverse momentum of the population of 
zero modes is of the same order as the Fermi momentum of the quarks. Eq. (77) is obtained 
by, as in Ref. [33], assuming that scattering involving the zero modes has the same relaxation 
time as quark-quark scattering in a non-superconducting medium (i.e. as if the flux tube core 
were infinitely large). We can then use the continuum quark-quark momentum relaxation 
time Ts from gluon exchange in a cold quark-gluon plasma (Ref. [32], eqn (28)) as a crude 
estimate of the relaxation time Tc for momentum transfer between bulk gapless quarks and 
zero modes inside the flux tube. 

Comparing (74) with (65) and (36) we see that the ratio of the zero mode force to the 
mutual friction force is /zm//mf ~ (i^oTcTr^a^)"^, assuming ~ 1. Using the estimates given 
above, 

. 0.003 i^^Y' C-^)' ( ^ f . (78) 

f^i V400 MeV/ V A / VO.Ol MeV/ ^ ' 

We conclude that the zero mode force is likely to be negligible relative to mutual friction. 



G. Timescale for expulsion of flux 

We can now estimate the time scale for the expulsion of the X magnetic field flux tubes 
from the 2SC core. As we noted above, there will be an outward force on the flux tubes 
at the point where they reach the nuclear mantle. The maximum force per unit length is 
given by (68), in which the energy costs of the magnetic field in the nuclear mantle have 
been neglected. The rate of outward movement of the flux tubes is given by balancing 
that force against frictional or pinning forces. There may be such forces arising from the 
nuclear matter, but we ignore them and only include the Aharonov-Bohm (mutual friction 
and lordanskii) and Magnus-Lorentz forces in the quark matter. Using (60), (65), (73), (72), 
we can see that the steady-state value of the vortex velocity is given by the force balance 
equation, 

pMi.{.Vs -Vl) xh + D\vl -vn) xn + r]{vL - v^) + /bf (r) = , (79) 

where /bf(^) is given by (68). We work in the reference frame that is uniformly rotating with 
the normal component (blue quarks and electrons) and we neglect possible small differential 
rotation between the superfluid and the normal fluid, so Vj^ = vs = in this frame. The 
lordanskii and Magnus-Lorentz forces then add to give a single transverse force 

-pVL X h + r]VL + fu{r) =0 , (80) 
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where p = puh — D'. From (72) and Sec. VIE we see that pml ~ P^ and D' ~ ap^, so we 
can neglect the lordanskii force and assume p ~ pml- 

We take the flux tube to lie in the z direction, and we calculate its position in the x, y 
plane using polar co-ordinates {r,6). We want to find r, the rate at which the flux tube 
moves outward. Solving (80) for the steady-state velocities f and 6, we find 

^ - Mr) , 



. ''V'^ (81) 

where fr is the radial component of the boundary force. We note in passing that f shows 
a non-monotonic dependence on the friction coefficient t]. As rj tends to zero one might 
expect the expulsion time to also tend to zero, and in the absence the Magnus-Lorentz 
force [p = 0) this would indeed be the case. However, in the presence of a non-zero Magnus- 
Lorentz force, the flux tube moves in an orbit around the center of the star, with the radially 
outward boundary force balanced by the resultant radially inward Magnus-Lorentz force. 

If the flux tube starts at radius tq at time t = and leaves the core (r reaches R) at time 
t = ti, then by solving (81) we find 



T 



21n(-) +1- 
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r/^ + p 
2ex V 
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(82) 



The factor in square brackets of order 1 for initial radii tq not too close to or R, so the 

v2,,3 



flux expulsion time for a typical flux tube is of order r. From (66) and (72), rj ~ « /ij, and 

„2,, d2 



p ~ yU^. So p ^ r], and using (14), (68), (36) we find 

Sa^^PgR^ 



T 



Tia^ In Kx 

Taking ~ 1, 



(83) 



The timescale for X flux tubes to be expelled from the 2SC core is therefore in the range of 
lO^'' years. 



VII. CONCLUSIONS 

Quark matter in the 2SC (or CFL) color-superconducting phase is a superconductor 
with respect to a broken "X" generator that is mostly color with a small admixture of 
electromagnetism. We have confirmed previous calculations [7] showing that quark matter 
in the 2SC phase will be a type-II X-superconductor if the quark pairing gap is above a 
critical value which is well within the expected range (8). Although the ambient magnetic 
field in the core of a neutron star is below the lower critical field for the formation of 
Abrikosov flux tubes containing X- magnetic flux, we argue that, when the quark matter 
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cools into the 2SC phase, the process of domain formation and amalgamation is likely to 
leave some of the X flux trapped in the form of flux tubes. The exact configuration and 
density of such tubes depends on details of the dynamics of the phase transition, but the 
density could be within an order of magnitude of the density of conventional flux tubes in 
proton-superconducting nuclear matter (11). Our calculations apply to 2SC quark matter 
in the temperature range Tisc < T <^ T2SC where T2SC is the critical temperature for the 
formation of the 2SC condensate, expected to be of order lOMeV (10^^ K), and Tisc is the 
critical temperature for self pairing of the blue quarks, which could be as low as 1 eV (10^ K). 

The 2SC phase contains three species of gapless fermions: two quarks ("blue up" and 
"blue down") and the electron. These are expected to dominate its transport properties. 
We do not discuss strange quarks, but our analysis is also applicable to phases with strange 
quarks present, as long as their pairing pattern does not break the Q gauge symmetry. 
Muons may also be present, but, like strange quarks, their higher mass gives them a lower 
Fermi momentum so they make a subleading contribution to the phenomena discussed here. 
We have calculated the Aharonov-Bohm scattering cross-section of gapless fermions with the 
X flux tubes (15), (36), and the associated collision (or relaxation) rate (55). A comparison 
with the collision time for Coulomb quark-quark scattering indicates that only at very low 
temperatures (T < 10'^ K or 10 eV) will the flux tubes dominate over thermal relaxation via 
Coulomb scattering. However, we defer a detailed calculation of the transport properties, 
including Coulomb and X-boson-mediated interactions, to future work. 

Because the ambient magnetic field in a neutron star is below the lower critical field 
required to force X-fiux tubes into 2SC quark matter, the trapped flux tubes will feel a 
boundary force pulling them outwards. We calculated this force for the case where the 
energy of the magnetic field outside the core can be neglected relative to the energy of flux 
tube. This force will be balanced by the drag force ("mutual friction") on the moving flux 
tube due to its Aharonov-Bohm interaction with the thermal population of gapless quarks 
and electrons (66), and also by the Magnus-Lorentz force (72). On this basis, we estimate 
that the timescale for the expulsion of X flux tubes from a 2SC core (84), is of order 10^° 
years. 

The work described here offers many directions for future development. 

(1) To get a full picture of the transport properties of 2SC quark matter one must calcu- 
late the relaxation rates associated with processes that do not include flux tubes, such as 
Q-Coulomb and X-boson-mediated interactions between gapless fermions. 

(2) We studied the regime where the cyclotron frequency is smaller than the thermal re- 
laxation time of the unpaired quarks (see Sec. VIA). It would be valuable to extend our 
analysis to higher magnetic fields and/or lower temperatures where the cyclotron frequency 
cannot be neglected. 

(3) It is important to resolve the disagreement in the literature over whether the Magnus- 
Lorentz force on flux tubes is cancelled by forces arising from the neutralizing background 
(see Sec. VI D). This is necessary for understanding flux expulsion from superconducting 
nuclear matter as well as more exotic flux tubes such as the ones that we described here. 

(4) We assumed that the X-flux tubes are stable, or at least have a lifetime that is long 
enough for them to play a role in transport. However there is no topological guarantee of 
their stability, and it is necessary to perform a calculation of their energetics, analogous to 
that of [13], and to investigate bound states on the string, which if present can enhance their 
stability [14]. 

(5) We focussed on the 2SC phase, but other phases may support flux tubes. The CFL 
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phase, which is the ground state of 3-flavor quark matter at asymptotically high densities, 
also has a gauge symmetry breaking pattern which resolves an external magnetic field in 
to an unbroken Q part, and a broken X part which could be carried in fiux tubes [46]. In 
this case also there is no topological guarantee of stabihty, and an analysis of the energetic 
stability is required. The CFL phase also features semi-superfluid vortices with non-zero 
magnetization [7, 47-50]. Since the CFL phase has no gapless charged excitations the asso- 
ciated phenomenology is likely to be quite different. In the CFL-KO phase there are charged 
kaon modes that can have an energy gap well below the pairing gap, so, if they have non-zero 
Aharonov-Bohm /3 factors, their scattering off flux tubes might be important. 

(6) We treated the thickness of the flux tubes as negligible, so scattering off them is domi- 
nated by the Aharonov-Bohm effect. In fact the thickness of the flux tube is comparable to 
the inverse Fermi momentum of the quarks (see (12)) and there will be flnite-size corrections 
to our results. Calculating them would require explicit construction of the radial proflle of 
the flux tube. 

(7) Some quark matter phases break the Q gauge symmetry. These include the 2SC phase 
at T < Tisc, and many other phases such as the color-spin-locked phase [16, 51]. It is 
interesting to ask what happens to magnetic flux in such cases: is the Q-superconductivity 
always type-I? (One suspects it may be because the gaps are usually small.) Will the dy- 
namics of the phase transition lead to trapped normal regions, and what is the tirnescale for 
their expulsion from the star? Could these phases retain X-flux tubes even after Q flux has 
been expelled? If X-flux tubes existed in a CFL core, for example they might experience 
the same sort of entanglement with superfluid vortices as is predicted in nuclear matter. 

(8) Neutron stars probably have layers of different phases. For a proper treatment of the 
dynamics of magnetic flux one would have to analyse how magnetic flux was connected 
between layers and pinned within layers, and the consequent additional forces on the color 
magnetic flux tube in a 2SC core. For instance, in a conducting nuclear mantle there would 
be eddy-current pinning of the magnetic flux; in a type-II superconducting and superfluid 
mantle there would be entanglement of nuclear Abrikosov flux tubes with superfluid vor- 
tices; and so on. There is also the possibility of different quark matter phases, such as an 
inner CFL core, inside the 2SC region. If it turned out that additional forces arising from 
these other regions of the star acted so as to allow expulsion of the flux tubes on a shorter 
timescale then this would have interesting astrophysical ramiflcations, such as a change of 
the magnetic moment of the star over this period of time. If the core contained a phase 
where X-flux tubes were entangled with superfluid vortices (as mentioned for the CFL phase 
above) then the rotational dynamics could also be affected. Observationally, this could pro- 
vide a new mechanism for glitches in neutron stars, since vortex-interface pinning force, 
derived above, may prevent a continuous flow of rotational vortices in the superfluid phases, 
in a manner analogous to vortex pinning in the crust [52] and the hadronic core-solid crust 
interface [53]. Other dynamical manifestations, such as, for example, the recently studied 
shear modes [54, 55] in the superfluid core and the post-jump relaxations (see Ref. [22] and 
references therein) will be affected as well. 
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Appendix A: Relating the scattering amplitude to Aharonov-Bohm cross-section 

The differential scattering probability is given by [56] 

dW = 2nS{e' - e)2n6ip. - P.)\MM'^ H (^^ffcr- (Al) 

It is assumed that a particle scatters off a heavy center; the momentum conserving deha 
function reflects the fact that there is no scattering along the vortex {i.e. , in the z direction). 
The quantities referring to the final states are primed, those to the initial state are unprimed. 
The a products is over all the final state particle phase space. The differential scattering 
cross section is 

dW . V eV 
da=—, J = - = — , A2 

3 v \p\ 

here j is the current, v the velocity, and V the volume. Substituting (A2) in Eq. (Al) we 
obtain for the scattering cross-section 

da = r'2n6{e' - e)2n6{p^ - p'.)\MM'^ U j^^^- (A3) 

The differential scattering cross section is obtained on writing d^p' = dp'j_p'j_d(j)'dz' (we 
restrict in the following the a product to one particle, since we consider elastic scattering 
and there is the same blue quark or electron in the final state). Thus, 

^ ' ^''^ ^ ' 2eV (27r)32£' 

The final state energy is e' = \Jp'1 + p'^, therefore e'de' = p'^dpl^. After integrating by 
means of delta-functions we obtain 

^ = j-i|M.,|2^— . (A4) 

The current density is given by j = v±/L'^, where f_L = [56]. Comparing (A4) with 

Eq. (15) we obtain 

\M,? ^ '-^^^ (A5) 
sm (0/2) 

Finally, substituting this result in Eq. (Al) we arrive at Eq. (53) of the main text. 



Appendix B: Phase space integrals 

Upon substituting the transition probability in the rate Eq. (49) and introducing mo- 
mentum transfer k = p' — p we find 

.0 = f ^«in^(^^) /(0[^-^]'^(^^~) (Bl) 
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where we replaced N^L/V = riy, which is the density of flux tubes per unit area. The 
integral Ip is defined as 



^JM [1 - Mp)\ mP) - eiU, - P))-^i^. 



(B2) 



The form of the scattering probability suggest that the phase space integral over d?p is 
convenient to carry out in the cylindrical coordinates by writing d^p = p±dp±d(f)dpz and 



= / dp±P±d4> . J' / dpjoip) [1 - /o(p)] -^S{e{p) - e{p - k^)). 

^[ZTiy J sin^(0/2) J e[py 



To do the inner integral note that in the low-temperature limit 

/o(p) [1 - /o(p)] = ^ T^4L= 6{p. -pI) 



de{p) 



pI-pI 



(B3) 



(B4) 



where in the last step we used cylindrical polar coordinates with p° = \/e% — p\. Carrying 
out the p^-integration we obtain 



T 



dp±p± 



4(27r)3 J p^^p^l 



where 



sin^(0/2) 



5{eF- \ pl + kl-2p^-kj_). 



(B5) 



(B6) 



We next specify the geometry of the scattering, by assuming that the vortex is along z-axis, 
vector p is directed along the x axis, and the scattering is in the x-y plane. If we denote 
the angle formed by the vectors p± and fc_L by x then, cosx = sin 0/2. We next note that 
identically dcj) = — 2(isin(0/2)/ cos(0/2) and define 



sm 



2pj 



(B7) 



The integral becomes 



Pf dsin{(j)/2) 



P±k± 
IGpfpI 



1 cos(0/2) [sin2(0/2)_ 

1 „A ki 



1 



5(sin 0/2 — sin 4>o/2) 



-.9 1 



2p± 



(B8) 



The integration limits were chosen for convenience [—1; 1]; the integral is then multiplied by 
a factor 2 to account for full 360 degree angle range. Inserting these results in Eq. (B5) we 
obtain 



2T 



dp±p] 



1 



1 



2pj 



pI-pD- 



(B9) 
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The integration is carried out using 



to obtain the final expression for the momentum integral 
Next we substitute this result in Eq. (Bl) and obtain 

The angular integral is equal vr, the remaining integral is (8/3)p|n, and we arrive at Eq. (54) 
of the main text. 
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